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q ' Abstract 

^ , Solutions of the general cubic complex Ginzburg-Landau equation comprising 

multiple spiral waves are considered. For parameters close to the vortex limit, and 
for a system of spiral waves with well-separated centres, laws of motion of the centres 
are found which vary depending on the order of magnitude of the separation of the 
centres. In particular, the direction of the interaction changes from along the line of 
centres to perpendicular to the line of centres as the separation increases, with the 
strength of the interaction algebraic at small separations and exponentially small at 
\ large separations. The corresponding asymptotic wavenumber and frequency are 

determined. These depend on the positions of the centres of the spirals, and so 
evolve slowly as the spirals move. 

1 Introduction 

> 

^ ■ The complex Ginzburg-Landau equation is one of the most-studied nonlinear models in 
physics. It describes on a qualitative level, and in many important cases on a quanti- 
tative level, a great number of phenomena, from nonlinear waves to second-order phase 
transitions, including superconductivity, superfluidity, Bose-Einstein condensation, liquid 
crystals, and string theory [3J. 

The equation arises as the amplitude equation in the vicinity of a Hopf bifurcation in 
spatially extended systems, and is therefore generic for active media displaying wave pat- 
terns. The simplest examples of such media are chemical oscillations such as the famous 
Belousov-Zhabotinsky reaction. More complex examples include thermal convection of 
binary fluids [13] and transverse patterns of high intensity light [S] . 
The general cubic complex Ginzburg-Landau equation is given by 

= y-h + ia) |^| 2 ^ + (1 + z6)V 2 ^, (1) 

at 

where a and b are real parameters and the complex field \& represents the amplitude and 
phase of the modulations of the oscillatory pattern. 

Of particular interest are "defect" solutions of (CD). These are topologically stable 
solutions in which \I/ has a single zero, around which the phase of \& varies by a non-zero 
integer multiple of 2n. When a = b these solutions are known as "vortices", and the 
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constant phase lines are rays emanating from the zero. When a 7^ b the defect solutions 
are known as "spirals" , with the constant phase lines behaving as rotating Archimedean 
spirals except in the immediate vicinity of the core. 

It is often convenient to factor out the rotation of the spiral, by writing 



This gives, on dropping the primes, 

di 
~dt 



i 6)^ = V 2 ^ + (1 - |^| 2 )^ + z#(l - k 2 - |^| 2 ), (2) 



where 



2n & 



1 + ba yv ' l + 6w' 
rotating single spiral waves are now stationary solutions of (J2j). The constant k is known as 
the asymptotic wavenumber, since it is easily shown that at infinity arg(^) ~ n<p±kr. An 
important property of spiral wave solutions is that the asymptotic wavenumber k is not a 
free parameter, but is uniquely determined by q [7] . Physical systems corresponding to (J2J) 
generally contain not one but many defects. Such complex patterns may be understood in 
terms of the position of these defects. Thus if the motion of the defects can be determined, 
much of the dynamics of the solution can be understood. 

Defect solutions of ([2]) behave very differently depending on whether q = (corre- 
sponding to a = b) or q 7^ (corresponding to a 7^ b). When q = the wavenumber 
k = 0, and a great amount is known about the solutions to (j2J). In particular, in a seminal 
work, Neu [9] analysed a system of n vortices asymptotically in the limit in which the 
separation of vortices is much greater than the core radius, using the theory of matched 
asymptotic expansions. By approximating the solution using near-field or "inner" ex- 
pansions in the vicinity of each vortex core and matching these to a far-field or "outer" 
expansion away from vortex cores, Neu derived a law of motion for each vortex in terms 
of the positions of the others, thus reducing to the solution of 2n ordinary differential 
equations (for the x- and ^-coordinates of each vortex). The interaction between defects 
in this case is long-range, essentially decaying like r -1 for large r. Neu's analysis has 
become the template for the analysis of the motion of a system of singularities in many 
equations, including more detailed models of superconductivity [10], [I], [5]. As we shall 
show, the key property of (j2J) that facilitated Neu's analysis is the linearity of the far- 
field equations. Thus the contribution from many vortices in the far field can be obtained 
by a simple linear superposition, and the motion of vortices is determined through the 
interaction of this far-field with the individual core solutions. 

When q > the wavenumber k > 0, and the situation is much more complicated, even 
for a single defect. Hagan [7] studied single spiral wave solutions of (J2J) and demonstrated 
that the asymptotic wavenumber k is uniquely determined by q. Using perturbation 
techniques he found asymptotic expressions for ip and the asymptotic wavenumber k as 
a function of q and the winding number of the spiral n. For small values of q a single 
defect has a multi-layer structure, with the solution comprising inner, outer and far-field 
regions in which different approximations hold. The transition from the outer region to 
the far-field occurs exponentially far (in q) from the centre (at what we shall term the 
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outer core radius), and the asymptotic wavenumber k is correspondingly exponentially 
small in q. This outer core radius is the radius at which the level phase lines switch 
from being essentially radial to essentially azimuthal. Thus for non-zero q there is a new 
lengthscale in the problem, with each spiral core having two lengthscales. In studying the 
motion of spirals it is no longer enough to say that they are well-separated by comparison 
to the core radius; now it must be determined whether the separation is large compared 
to the inner core, but small compared to the outer, or whether the separation is large 
compared to the outer core radius so that the interaction is truly far-field. When the 
separation lies between the inner and outer core radii the interaction is algebraic, but 
when the separation is large compared to the outer core radius the interaction of defects 
decays exponentially 

The fact that the outer equation for the phase of ip is nonlinear, so that the con- 
tributions from multiple defects may not simply be added, along with the exponential 
scaling of the outer variable, explains the difficulty in applying Neu's techniques to the 
general case of non-zero q. Thus, despite much work and some partial results [21 H2], the 
interaction of defects in the case of non-zero q was not completely understood. However, 
recently in [1] a set of laws of motion for N spirals with unit winding number was derived 
systematically in the limit 0<g<l. The aim of the present work is to give the details 
of that calculation. 

We start in Section [21 with reviewing the general asymptotic scheme that determines 
the spirals' mobility when q = 0. In Section [3] we inspect the equilibrium solutions to 
(J2J) for q > formed by a single spiral, highlighting the existence of two distinguished 
outer regions where either the azimuthal or the radial components of the phase function 
dominate. These two sections serve as a template for the analysis of multiple-spirals 
patterns, which is performed in Sections H] and El In Section H] we derive a law of motion 
for spirals which are separated by distances comparable to the outer core radius; since this 
is a distinguished limit we refer to it as the canonical separation. The interaction at the 
canonical separation is found to be exponentially small and it takes place in the direction 
perpendicular to the line of centres of the spirals. In Section [5] spirals are assumed to be 
separated by distances lying between the inner and outer core radii; we denote this as the 
near-field separation. In this case the interaction becomes algebraic with a component 
along the line of the centres. 

Changing the separation of the spirals while keeping the parameter q fixed is equivalent 
to varying q at fixed separation. We will see that the near-field separation is needed to 
interpolate between the canonical separation and the case q — 0. 

2 Interaction of vortices in the Ginzburg-Landau equa- 
tion with real coefficients 

Without lost of generality and to simplify the calculations we consider equation ([2D with 
b = 0. We return briefly to the case of general b in the appendix. For q = this reads 




(3) 
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which by writing \I> = fe %x with / and x rea l an d separating real and imaginary parts we 
may write as 

/* = V 2 /-/|Vx| 2 + (l-/ 2 )/, (4) 
f X t = 2V X V/ + /V 2 X . (5) 

We wish to determine the law of motion for well separated vortices following [9]. We 
assume that the separation is 0(e _1 ), with eCl. This leads to an "outer region", in 
which x is scaled with e _1 , and an "inner region" in the vicinity of each vortex. Matching 
the asymptotic expansions of the solutions in each of these regions leads to the law of 
motion for each vortex centre. 



2.1 Outer region 

In the outer region we rescale x and t by setting X = ex and T = e 2 fit; here fi is a small 
parameter (the timescale for vortex motion) which will be determined later. We will find 
that \x is logarithmic in e. With this rescaling (jlj), §5§ read 

eV/r = e 2 (V 2 /-/|V X | 2 ) + (W 2 )/, (6) 
fifXT = 2Vx-V/ + /V 2 X . (7) 



Expanding in powers of e as 



/ ~ fo + e 2 h + 
X ~ Xo + e 2 Xi + 



we find 



fo = 1, (8) 

®Xo ^2 f n \ 

V-Qf = V xo (9) 



Now expanding xo f° r small \i as 



Xo ~ Xoo + fiXoi H 

and substituting into @ we obtain at leading order 

V 2 Xoo = 0, 

with solution 

N 

Xoo = X! n ^' 

i=i 

where <fij is the polar angle measured from the centre of the j th vortex Xj, and nj is the 
winding number (or degree) of the j th vortex. At the next order in \x we find 

V Xoi - -jyjT, 
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with solution 

N 



2 ^-j~j~ ~j~ v j dT 



where Rj is the distance from the jth vortex. Continuing to 0(/i 2 ) we find 

1 N ( 

.7=1 V 



(DC A / rfX 7 

'Or ) v rr Hr 



In general we find that Xom is 0(R m log J?) as i? — > 0. 
2.2 Inner Region 

We rescale near the £th vortex by setting X = + ex to give 

e ^( e/T "^" v/ ) = v2 /-/l v ^l 2 + ( 1 -/ 2 )/, (10) 

.cDQ \ 2 



^UfxT-f-^- -Vx) = 2V X V/ + /V i x, 

or equivalently 

e/ife^T-^-V^) =^(1-|^| 2 )+V 2 ^. (12) 



dT 

Expanding in powers of e as \& ~ \l/o + + • • • we find at leading order 

V 2 * + *o(l-|*o| 2 ) = 0. (13) 
This is just the equation for a single static vortex, with solution 

*o = /oe !M = /o(r)e'W +c ™, (14) 

where 

d 2 f ldf f n\ 3 
"77 H j Jo - 7 + Jo — Jo — u ) 

/(0) = 0, / - 1 as r - oo, 

and C(T) is determined by matching. It is well known that this equation has a unique 
increasing monotone solution [6]. 

Continuing with the expansion we find at first order in e that 

- /iW ■ = V 2 ^i + *i(l - |* | 2 ) - M*o*l + *o*0- (15) 
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2.2.1 Solvability Condition 

Since the linear operator on the right-hand side of equation (1151) is self-adjoint, and the 
homogeneous version of (IToT) is satisfied by the partial derivatives of \&o ( as can be seen 
by differentiating (TO]) ), by the Fredholm Alternative there is a solvability condition on 
(|15p which we can write as 

-Is JA*-*®- 1 *^} "' 

where d is an arbitrary constant vector and D is an arbitrary region in the plane. Taking 
D to be a ball of radius r we find, after some calculations, 

- of ■ d ) i ( s <« 2 + -4) * - - r & + t) * ■ d * < i6) 

where x ~ Xo + e Xi + ' ' " ■ After matching with the outer region to determine xi> equation 
f TT6l) will determine vortex velocity dX-g/dT. 



2.3 Asymptotic matching 
2.3.1 Inner limit of the outer 

We express the leading-order (in e) outer xo i n terms of the inner variable x given by 
X = + ex so that Ri = er, <fie = 0. Then, Taylor-expanding xo, 

Xo ~ Xoo + UXoi + • • • 

- G(X,) - n t j log er e^ • + eVG(X,) ■ x + 0(e 2 ), (17) 



where 



G(X) = E^-^Efl X - X ^ lo gl X - X il^ 



2.3.2 Outer limit of the inner 

The leading-order phase in the inner region is 

Xo = n e( p + C(T). 

This matches with the first term in ( fTTj) providing C(T) = G(X^). 

The two-term inner expansion for the phase is Xo + e Xi- This should match with the 
two-term inner expansion of the outer (j!7p . Since \x = 0(l/logl/e) and all logarithmic 
terms need to be matched at the same time, to perform this matching we need to take the 
full fi expansion of both terms. Fortunately only the expansion of G(X^) and VG(~Xg) 
involve infinitely many terms in /i, and these are evaluated at X^ and are therefore 
independent of x. We see that the expansions match if 

Xi n^logere^ ■ -Jj^// + VG(X^) x (18) 

as r — > oo, which is good to all orders in [i. 
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2.4 Law of motion 

We can now use the matching condition (118j) in the solvability condition (1161) to find a 
law of motion for each vortex. 

As r — > oo the left-hand side of (TT5]) is 

-/i7r(n^logr + a)d- (19) 



where 

a = lim 



£ (s(^) 2 + ^) ds-\ogr , 
is a constant independent of fi and e. Using (Tig]) in the right-hand side of (|TU|) gives 

- im{\j2 + \ogr)nl d • ^ - /m| logeTrd • ^ + 2vrd ■ VG ± (X,), (20) 

where - 1 represents rotation of a vector by ir/2. Since /i is a small constant, we find that 
the only way to make the left and right hand sides balance is by taking fi = 1/ log(l/e), 
giving 

-a,d-§ = ^.§-4d-^ + 2M-VGU). (21) 
Since d is arbitrary this can be rearranged to give the law of motion 

dX > 2n,VG-(X,) 



<ff n 2 + /i(a -nf/2) 



where 













[dT ^ 



+ 0(/i 2 ). 

Note that the expression ( 12~2j) is accurate to all orders in \x\ the expansion in /i is necessary 
only to evaluate VG. To leading order in /i the law of motion reads 

(23) 

dT m^lXe-XA' v ; 



In particular, for two vortices at positions (X\, 0) and (X 2 , 0) with X\ < X2 the laws 
of motion are given by 

dXl 2 .OVOGO. (24) 



c?T rii \Xi — X 2 

dX 2 m / 2 



,0 +0(Ai). (25) 



dT n 2 \X 2 - X 

The direction of motion is always along the line of centres, with like vortices repelling 
and opposites attracting. 
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2.5 An alternative matching procedure 

The analysis above is based follows the method that was used in [9]. However, when 
q > there is the added complication of three small parameters, /i, e and q, rather than 
just two, /i and e. This has several implications on the way we compute the asymptotic 
expansions, and affects the way the matching procedure must be carried out. 

In the analysis above the key step which determines the law of motion is matching 
the first-order inner solution xi with the correction to the outer phase due to the other 
vortices. When q = we are fortunate that we do not need to expand the inner xi i n 
powers of fi, so that we in effect retain all orders of [i while matching in e. For q > we 
are not so lucky, and we need to expand both the inner and outer solutions in powers of ji 
to make any progress. Since we then no longer have the full n expansion of either region, 
we cannot match them together. We can get around this deficiency by writing down and 
solving equations for the outer limit of the inner solution. By solving these equations we 
are able to resum the infinite series in fi that is present in the inner region, and this is 
exactly what we need when matching. To illustrate this new procedure we apply it here 
to the q = case, where we check the results against the known solution above. 

2.5.1 Outer limit of the leading-order inner 

Rather than solving the inner equations and writing the solution in terms of the outer 
variable and expanding, we rewriting the leading-order inner equations in terms of the 
outer variable R = er to obtain 

= e 2 (V 2 /o-/o|V X o| 2 ) + (l-/o 2 )/o, (26) 
= e 2 V-(/ 2 V X o). (27) 

We now expand in powers of e to give the outer limit of the leading-order inner solution 

as 

Xo ~ Xoo + e 2 Xoi H 

fo ~ Zoo + ejoi H • 

Substituting these expansions into ( 1261) . ( 1271) gives 

Too = 1, /oi = -^|Vxoo| 2 , (28) 
= V 2 £oo, (29) 

with solution xoo = nk4>. 

2.5.2 Outer limit of the first-order inner 

We write the first-order inner equation in terms of the outer variable to give 

-e/i^-V/o = e 2 V 2 / 1 -e 2 / 1 |Vxor-2e 2 / Vxo- Vxi + Zi-3/ 2 /!, (30) 
-fief^-Vxo = 6 2 V-(/ 2 Vxi) + 6 2 V-(2/ /iVxo). (31) 
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We now expand in powers of e as 



Xw . - 

Xi ~ + 

e 

/i ~ /io + e/n + 



to give 



/io = 0, 

/n = -Vxoo • Vxio, 

V72- V7^ 1 /X . 

V Xw = -^—^7 ■ VXoo = -n k n—j^- ■ — = n fc — (Ki sm0 - K 2 cos. 



where we have written 



Thus 

fiRXogR fiR\ogRdX k 
Xw = n k (Kisin0- K 2 cos0) = -n k — ■ 

plus a homogeneous solution which comes from matching with the outer. We see that 
this homogeneous solution is 

X- VG(X fc ), 

giving 

fiRlogR dX k V7^rv n ^qo^ 

Xw = ~n k — -e + X- VG(X fc ). (32) 

This expression is the outer limit of the full \i expansion of the first-order inner solution. 
Rewriting in terms of the inner variable gives 

fir log edX k vt^/y n , 

Xi n k — — ■ e + x ■ VG(X fc ) + • ■ ■ , 

which is ( I18p as expected. Note that first term in (1321) . which is the particular integral, 
was obtained previously from the O(fi) terms in the outer solution; this time we have 
used only the leading-order outer solution. Thus this method allows us to match with the 
inner solution when we only know a few terms in the /i-expansion of the outer solution. 



3 Equilibrium spiral wave solutions 

We now consider equation (J2J) and analyse equilibrium spiral wave solutions in the limit 
where the parameter q is small. These equilibria, which correspond to single spirals with 
arbitrary winding numbers, were studied by Hagan [7] who showed that the asymptotic 
wavenumber k, and thus the frequency of the corresponding periodic solution to (JT]), 
is uniquely determined by the parameter q. We shall recast Hagan's results in a more 
systematic asymptotic framework in a way which generalises to the many-spiral solutions 
considered in SfH an d ill 

With ip = je %x and / and x real, we seek solutions of the form / = f(r) and x — 
n<f) + ip(r) with f(r) and <p r { r ) bounded as r — > oo and 

^(0) = 0, /(r) ~ Cr n as r -> 0, 
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where C is some positive constant, and a subscript denotes partial differentiation. In fact 
it can be shown [7j that bounded solutions of (J2J) satisfy f(r) — > (1 — A; 2 ) 1 / 2 and (p r — > —A; 
as r — > oo. 

We start by introducing an auxiliary parameter e. For the multi-spiral case e will 
represent the inverse of the spiral separation; here it represents the inverse of the outer 
core radius, which will be defined shortly. Rescaling (j2J) onto this new (outer) lengthscale 
by setting X = ex gives 

Z6 2 Q? 2 

= e 2 V 2 ^ + (l+zg)(l-H 2 )^ i), (33) 

q 

where we have introduced the new parameter a = qk/e. This can now be seen as an 
eigenvalue problem for a(q) which provides the relationship between k and q. The outer 
core radius of the spiral is the value of e which makes a of order one as q — > 0. 



3.1 Outer Region 

With / = f(r) and \ — n 4> + ^( r ) equation ( 1331) becomes 

= e2 (/" + ^)- e2 /(^ + (^) 2 )+(l-/ 2 )/, (34) 

o = f(/W)' + g/ 2 (W 2 )-— / 2 - (35) 

it q 
Expanding in powers of e as 

/ ~ / (X;g) + e 2 / 1 (X;g) + --- , 
if ~ Vo(X; q) + e 2 ifi(K; q) H , 

we find 

/o = l, h = -\ (jj^ + (v') 2 ) • 

Equation flHBj) is a Riccati equation and can be linearised through the transformation 
<p = (l/q) logFo to give 

^ + f + ^o(^-« 2 )=0, (37) 

with the general solution H (R) = K inq (aR) + XIi nq (aR) where A is an arbitrary real 
number, and K inq and Ii nq are the modified Bessel functions of the first and second kind. 
Only when A = is the function ip monotone [7], so that the iso-phase contours are 
spirals. Then 

Xo ~ ncf) + - \og{K inq (aR)). (38) 
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3.2 Inner region 

We return to the inner (original) scaling by setting X = ex to give 

o = v 2 /-/|Vxl 2 + (i-/ 2 )/, 

= V-(/ 2 Vx)+g(l-/ 2 )/ 2 -^. 

Expanding 

/ ~ fo + e/i + e 2 f 2 + ■ ■ ■ , 
Lp ~ (p + e^i + e V2 H 

the leading-order equations are 

fS+& - /o + K) 2 ) + (1 - /oVo = 0, 

/o(^o + y)+2/^ + g(l-/ 2 )/o = 0. 
We now expand the leading-order solution in e in powers of the small parameter 

fo ~ foo + fo2Q 2 + /o4<? 4 H , 

V^oo s 
<£>o ~ h <A)2<7 + </2o4<? H • 

Substituting these expansions into ( 1391) and ( 1401) and equating powers of q gives 

V9 o = -Do, 
= ZoU^-- 2 ^ + (Wo 2 o)/oo, 



i /" r 

^02 = -772-/ s/ooC 1 - /oo) ds, 
'ioo </0 



'00 ^0 

with boundary conditions 

/oo(0) =0, lim /oo(r) = 1, 

r— >oo 

where -Do is a rea l constant to be determined by matching. 

3.3 Asymptotic matching 
3.3.1 Outer limit of the inner 

From expressions (1441) - (1451) we find that as r — > oo 

, 2 logr + c n 
^02 ~ ~n + ■ • • , 

where constant given by 



lim 



\ 2 ([ foo(s) (l - foo(s) 2 ) sds-n 2 log(r)) . 
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However, since we will see that q is logarithmic in e, in order to match with the outer 
expansion we need the outer limit of the full expansion in q of the leading-order inner 
solution in e. This is found to be of the form 

1 oo 

/o ~ 1 " ~2 E «;{^ 2 0og(r) + c n )} 2 * + • • • , (48) 

i=0 

oo 

^ ~ -- $>^>g( r ) + c «)> 2t+1 + • • • > ( 49 ) 

i=0 

where a, > and /3j > are constants independent of q and n. The necessity of taking all 
the terms in q when matching can be seen, since expansions (I48p and (l4T?j) are asymptotic 
provided g(log(r) + c„) <C 1; with g of order l/log(l/e) and r of order 1/e all the terms 
are the same order. 

In principal, when matching logarithmic expansions, the full logarithmic expansion 
in both the inner and outer regions needs to be found. However, when these are finally 
compared term by term they must both be written in either the inner or outer variable. 
Thus, in fact, it is enough to have the full logarithmic expansion in one of the two regions 
only, providing the expansions are written in terms of the other variable when comparing 
terms. 

For the single spiral case we have the full expansion in q of the leading-order outer 
solution (1381) . but only the first few terms in the g-expansion of the leading-order inner 
solution (T4*3l) -( r45l) . Thus we can complete the matching by writing in outer solution ( l38l) 
in terms of the inner variable, expanding in e and q, and comparing terms with ( I43l - fj45l) . 
However, for the multiple-spiral case we will not be able to solve to outer problem to all 
orders in q, so such an approach will not be available. 

Fortunately there is a method which allows us to sum the q expansion of the inner 
solution (|43 l) -( r45l) . and this method is still applicable in the multispiral case. As in §2.51 
rather than solving the leading-order inner equations, writing the solution in terms of 
the outer variable, and then expanding again ready to match, the trick is to write the 
equation for the leading-order inner solution in terms of the outer variable, and only then 
solve it. Hagan [7] thinks of this as a middle region expansion, but actually we are just 
writing down the equations satisfied by the outer limit of the inner expansion. 

Thus we begin by rewriting the leading-order inner equations ( 1391) and (T40l) in terms 
of the outer variable R = er to obtain 

= e 2 (^ + |-/o(^ + K) 2 ))+(l-/o 2 )/o, (50) 

= ^(RfM' + Qfo^-fo)- (51) 

Note that these are similar to, but not identical to, the equations in the outer region 
(|34l)-(l35l). We now re-expand in powers of e as 

¥o ~ <Poo(q) + e 2< ?oi(<?) H , (52) 

fo ~ /oo(<?) + e 2 /oi (<?) + ••• • (53) 

The leading-order term in this expansion, tfoo{q), is just the first term (in e) in the 
outer expansion of the leading-order inner solution, but now it includes the resummed 
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contribution from all the terms in q. Substituting ( !52l) - (j53l) into (l50l) - (l5"Tl) gives 



2 \ 2 



Zoo - 1, foi - yjp + (^oo) 2 
^o + ^ + g(^ + (^oo) 2 )=0. 



(54) 



Equation (1541) is again a Riccati equation which can be linearised with the change of 
variable (poo = (V?) logS to give 

* S+ R +i W So = > (55) 

with solution 

H = A(q)e- iqn R iqn + B{q)e tqn R~ iqn , (56) 

where A and S are constants that depend on q and the factors e ±tqn have been included 
to facilitate their determination by comparison with the solution in the inner variable. 
One relationship between A and B is determined by writing </? o in terms of r, expanding 
in powers of q, and comparing with 049p . A second is given by matching with the outer 
solution. Expanding the constants in powers of q as 

A(q) ~ -Aq + A 1 + qA 2 + ■ ■ ■ , (57) 
Q 

B(q) ~ -^S + B 1 + qB 2 + ■ ■ ■ , (58) 

writing Hq in terms of r, and expanding for small q we find 

% = A(q)e iqnlosr + B(q)e- iqnlosr 

1 „ , , \ A . , g 2 n 2 log 2 r 

A + A\ + qA 2 H 1 + zgnlogr 



+ -S + B x + qB 2 + ■ • ■ (l- iqnlogr- - —2— + 



Aq + Bo 



A t + B 1 + (A - B )in\ogr 



q 

+ q[A 2 + B 2 + {A x - B^in log r - n 2 ^— ± ^ log 2 r 



so that 



H' {r) (A -B )ni ({Ai-B x )in 2 \ogr 



^° (r) qHo(r) ~ r(A + S ) + ^ V (A + S )r 



2 (A -S ) 2 logr _ ^(A -S )(A 1 + S 1 ) n 
n (A + S ) 2 r (Ao + S ) 2 
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Comparing with ( 1431) . and (I4T1) we see that 

A) - B =0, (59) 

1 T 5 z = - nc «' ( 60 ) 

The remaining equations determining A and S will come from matching with the outer 
region. 

3.3.2 Inner limit of the outer 

Since we have now been able to sum the g-expansion in the inner, we do not need the 
full g-expansion of the outer, but only the leading-order term. Since we have used the 
same transformation from to H in the inner and the outer, it is easiest to perform the 
matching in terms of H . We find 

H {R) = K mq {aR) ~ K (aR) + 0(q 2 ). (61) 

Thus the inner limit of the outer is 

H (R)~ -log — - 1 + 0(q 2 ) (62) 

where 7 is Euler's constant. This should match with H (R) given by (|5*B1) . Exanding 
H (R) in powers of q gives 

Aoe -i<,nlo g e + Boe i q nlo g e 

H ~ 1 (63) 

q 

Comparing flgH} and (E2D and using fl52J we find that e iqnloge + e - iqnl ° se = (D(q), so that 

g|nloge| = — + vq + 0(q 2 ). (64) 
Under this condition (l6"3"j) becomes 

^0 ~ -\n\(A + Bo) \ogR- (A Q + B )v + zsign^)^ - 
Comparing with ( 1621 we see 



H(A + B ) = 1, 
zsign(n)(Ai - B x ) - (A + 5 )v = - log a + log 2 - 7. 



Eliminating A and I? using (jo9j) and (160]) gives, finally 

a = 2e c "+^M-T, ( 65 ) 

so that the eigenvalue a is now determined, and is indeed an order-one constant. Recalling 
that a = qk/e and using (|64|) we find that the asymptotic wavenumber is given by 

k(q) = - exp ( c n - 7 - j (1 + o(l)) (66) 
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in agreement with Hagan [7] . The corresponding frequency of the pattern is given by 

u = q(l- eV) = q- 4 qe **-*r-*Mn\ (i + (1)). ( 67 ) 

We see from the expression for the phase in the outer region that the outer core radius 
corresponds to the radius at which the iso-phase lines have become essentially azimuthal. 
We will see that for the many-spiral case this corresponds to the interaction between 
spirals becoming perpendicular to the line of centres. For separations less than the outer 
core radius the interaction is algebraic, but for larger separations the interaction rapidly 
becomes exponentially small. Thus the outer core radius forms an effective "region of 
influence" for each spiral. 

4 Interaction of spirals at the canonical separation 

We now want to combine the methods presented in $2] and Sj3] to consider the interaction 
of well-separated spirals when < g < 1. In each of $2] and $3] we had two parameters 
to relate to each other (e and /i, and e and q respectively); here we will have to consider 
the relative sizes of all three parameters e, \x and q. We saw in $3] that there is an outer 
core radius at which H switches from logarithmic growth to exponential decay. Since the 
interaction between spirals occurs as a result of the phase, it is clear that the interaction 
will crucially depend on the relative sizes of the spiral separation and the outer core 
radius. 

We start by considering the distinguished limit in which the spirals are separated by 
distances of the same order as the outer core radius as q 0. Since this outer core radius 
varies exponentially with winding number, we now make the assumption that all winding 
numbers are ±1, so that all spirals have the same outer core radius. This means that we 
again have a = kq/e of order one, and a spiral separation of order e K ^ 2q . 

We will show that the spirals interact, to leading order, in the direction perpendicular 
to the line of centres, with a velocity of order e 2 /| loge| (as in We will find that for a 
pair of spirals, in contrast to fluid vortices, the direction of motion of each depends only 
on its own winding number, and not on the winding number of the other spiral. 

4.1 Outer region 

Since we are considering spirals separated by distances large compared to the inner core 
radius, the solution in the vicinity of each spiral will be a small perturbation of the 
single spiral solution. The whole pattern will have an associated frequency or asymptotic 
wavenumber. However, in the multi-spiral case, as we shall show, this frequency is no 
longer constant but varies slowly as the spirals move, tending to the single spiral limit as 
the separation tends to infinity. Thus we need to remember that when we use equation 
(J2J) the eigenvalue k (or equivalently a) may depend on time. 

As before we rescale time and space in (J2J) by setting X = ex, T = fie 2 t, to give 

e 2 ^ T = e 2 V 2 iP + (1 + iq)(l - |^| 2 )^ - (68) 
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Writing = fe lx as usual and separating real and imaginary parts in (168j) gives 

/ie 2 /T = 6 2 V 2 / - 6 2 /|V X | 2 + (1 - f)f, (69) 
^Pxt = 6 2 V • (fV X ) + g/ 2 (l - f) ~ e 2 -f. (70) 

q 

Expanding in powers of e as 

/ ~ / (X,T;g, / i) + e 2 / 1 (X,T;g, /.) + ••■ , 
X ~ Xo(X,T;g,/i) + e 2 Xi(X,T;g, //) + ■■• , 



we find 

/o = l, /i = -^|V X0 | 2 , 

a 2 

VXot = V 2 Xo + g|Vxo| 2 • (71) 

9 

It is tempting now to linearise (1711) via the Cole-Hopf transformation xo — (1/?) log/i as 
in $3], to give 

M /i T = V 2 /i - ct 2 /i. (72) 

Then, since the equation is linear, we could sum up the contributions from each spiral to 
give 

JV 

h = ^(T)e^K iqn (aRj), (73) 

where Rj and <pj are the polar variables centred on the j th spiral and the weights /3j 
depend on the slow time variable T. This function has the right type of singularities to 
match with the spiral core when we expand it locally. Unfortunately, when the transor- 
mation is undone to return to Xo an d ip, we see that ip is no longer single-valued. This 
problem with multivaluedness is pointed out in [UJ and [12], where the authors rightly 
claim that it invalidates the use of the Cole-Hopf transformation. Nevertheless, in what 
follows we will show that this transformation can be used to advantage without causing 
ip to become multivalued, providing care is taken. The key is the observation that for a 
single spiral the dependence of x on occurs at 0(1), not 0(l/q), so that at leading order 
the Cole-Hopf transformation can be used without difficulty Then, at first order, the 
single-valueness of ip will be maintained by introducing exactly the right multivalueness 
in h. 

To simplify the exposition we now make the assumption, which will be justified a 
posteriori, that /j, and q are the same order as e — > (we will see that both are 0(1/ log(e)). 
We therefore write \i — qjl and treat ft as 0(1). Expanding xo m powers of q as 

Xoo . 

Xo ~ h Xoi H 

q 

gives, to leading order, 

= V 2 Xoo + |Vxoo| 2 -« 2 . (74) 
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Linearising (174j) through the Cole-Hopf transformation xoo = log ho gives 

= V 2 /i - a 2 h . (75) 

Now, at leading order, the solution for a single spiral is given by ( J6T1) ; for multiple spirals 
we can some these to give 

N N 

h = ^ (3j (T)K (aRj), xoo = log^/3,(T)iT (^). (76) 

3=1 3=1 

Note that, because the leading-order solution does not depend on 0, there is no problem 
with multivalueness of ip. The weights /3j will be determined by matching with the inner 
expansion. 



4.2 Inner Region 

We rescale near the I th vortex by setting X = + ex to give 

^(Vr-^-V/) = V 2 /-/|Vx| 2 + (l-/ 2 )/ (77) 



^f 2 [e X T-^-Vx) = v-(/ 2 v x ) + g(i-/ 2 )/ 2 -^ (n 



or equivalently 



Expanding 



dX \ e 2 a 2 
^t~^- = V 2 V + (1 + iq)(l - H> - i— */> (79) 



/ ~ / (x;g,/i) + e/i(x;g,/i)H , 

X ~ Xo(x;g,/i) + exi(x;g,/i) H , 

the leading-order equations are 

= V 2 / -/o|Vxo| 2 + (l-/o)/o, (80) 

= V-(/ 2 V X o)+g(Wo 2 )/o, (81) 

or equivalently 

= V 2 ^o + (1+ ^o(l -W 2 ). (82) 

The effect of the other spirals is felt at higher order, so that the leading-order inner 
solution is of the form / = fo(r) and xo = nefi + V 9 o( r '), where 



/o+^o-/o(^ + K) 2 )+(l-/ 2 )/o = 0, 



53) 



fo[^o + -) +2/^ + g(l-/ 2 )/ = 0, 



r 
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where are of course the equations for a single spiral (1391) - (1401) . Expanding in powers of q 
as in (I4"l]) - (l4"2]) we find that (p 00 — <Poo(T), with / o and <^ 02 satisfying (TJlj)-(j46]). 
At first order in e we find 

- ^ ■ Wo = VVi + (1 + - l^ol 2 ) - tf>ofrM* + rM) (85) 

or equivalently 

-/^■V/o = V 2 /i-/i|Vxo| 2 -2/ Vxo- Vxi + /i-3/ 2 /i, (86) 

-/i/ 2 ^-Vxo = V-(/ 2 Vxi) + V-(2/ /iVxo) + 2g/o/i-4g/ 3 /i. (87) 



The solvability condition on (|85j) will give the law of motion. We derive this condition 
in §4.41 First, we match the leading-order solutions in each region to determine the 
weights (3j and the eigenvalue a. 

4.3 Asymptotic matching 

4.3.1 Outer limit of the leading-order inner 

We use the same trick as in Section [3] to sum the g-expansion of the outer limit of the 
leading-order inner solution. We begin by rewriting the leading-order inner equations 
( I80p and (181H in terms of the outer variable Rg = er to obtain 



= e 2 (V 2 / -/o|V X o| 2 ) + (l-/o 2 )/o, (88) 
= e 2 V -(/ 2 Vxo) + g(l-/o 2 )/o- (89) 

We will now denote Re by R to simplify the notation. We now expand in powers of e as 

Xo ~ Xoo(g) + e 2 Xoi(<?) H , 

fo ~ /oo(g) + e 2 /oi(g) H — • 

The leading order term in this expansion Xooil) is just the first term (in e) in the outer 
expansion of the leading order inner solution, including all the terms in q. Substituting 
these expansions into (1581 . (159!) gives 



foo — 1, foi — — 2 I^Xool 2 , 
= V 2 Xoo + g|Vxoo| 2 . (90) 

Equation ( 1901) can be linearised with the usual change of variable xoo = (1/?) log/^o to 
give 

V 2 h = 0. 

As in the single spiral case the relevant solution is of the form £00 = n e4> + <p{R) so that 
h = e qni<f e w{R) = e qn "f>H (R) where 

^o' + f + g 2 § = 0, (91) 
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with solution 

H = A e (q, T)e~ iqne R iqne + B e (q, T)e iqn * R~ iqne , (92) 

where the constants An and Bg may depend not only on g but also on the slow time T, 
and may be different at each vortex; as before the factors e ±iqrLl have been included to 
facilitate comparison with the solution in the inner variable. Expanding 

A e (q) ~ — + An + qA i2 + • • • , (93) 
Q 

B,{q) ~ — + B ei + qB n + ■■■ , (94) 
Q 

writing x o i n terms of r, expanding in powers of q, and comparing with (l4"3~j) . and (|4"7|) 
as in the single spiral case gives 

A m - B m =0, (95) 
(A ei - B n ) 



Aiq + Bg 



-n e c ne for £ = 1,..., N. (96) 



The remaining equations determining Ag and Bg will be fixed when matching with the 
outer region. 

4.3.2 Inner limit of the outer 

To compute the inner limit of the leading-order outer solution we rewrite solution ( 1761) 
in terms of the inner variable by setting X = + ex and expand in powers of e to give 

N 

j=i 

P t log — - fa + G(Xt) + ex ■ VG(X £ ) + • • • , (97) 

where 

N 

G(X)= Yl ^( r )^o(«|X-X y |). (98) 

4.3.3 Leading order matching: determination of the frequency 

We now match the inner limit of the leading-order outer solution (pTI) with the outer limit 
of the leading-order inner solution (1921) . Since we have summed the g-expansion of the 
outer limit of the leading-order inner, while we have only the first term in the g-expansion 
of the inner limit of the leading-order outer, we must write both expansions in the outer 
variable R before comparing terms. 

The g-expansion of the inner limit of the leading-order outer solution is 

Xo „ l^o + . . _ I log fp tl ?R_ p + G{Xe) \ + . . . (99) 

q q V 2 J 
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The outer limit of the leading-order inner solution is 

Xoo = n £ (f)+ - log H (R) 
Q 

= n t </>+- log (A e (q)e~ iqne R iqne + B e {q)e iqne R~ iqne ) 

= ri£0+-log H . (100) 

9 V 9 / 

Comparing (11001) and (I9"9l we see that, as in the single spiral case, we require 

e iqn e lege + ^iqn.loge = (101) 

so that, since n e = ±1, 

q\ogl/e=^ + uq + 0(q 2 ). (102) 

This is the relationship between q and e required for a to be of order one, and is equivalent 
to assuming that the typical spiral separation 1/e = (^(e 71 "/ 29 ). 
Assuming (I102p holds, equation (11001) becomes 

£oo ~ - log (-(Aeo + Beo) log/2 - (A i0 + B m )v + i sign(n £ )(A £1 - B a )) + • ■ ■ . (103) 
9 



Comparing with ( 1991) gives 

+ = A, 

iBign(n/)(i4a-S/i)-(^o + S«))i/ = -/%(loga - log2 + 7) + G(X*). 
Eliminating and £^ using ( |95l) and ( |96l) and ( |98l) gives 

AT 

- ( Cl + = -A(log« - log2 + 7) +^2/3jK {a\Xt - Xj\). (104) 

Since ( I104p holds for each spiral this is a system of N linear equations for the unknown 
weights (3j, j = l,...,n. Since the system of equations is homogeneous, a non-zero 
solution will exist only if the determinant is zero: this is the condition which determines 
the eigenvalue a in the multispiral case. Note that the weights /3j and the eigenvalue a 
(and therefore k) depend on the position of the spiral centres, and will therefore evolve 
on the slow timescale T. 

For the case of just two spirals equations (I104f) become 

-{ci+u)f3 1 =-/3 1 (loga-log2 + 7 )+^o(a|X 1 -X 2 |), (105) 
-(ci + = - /3 2 (loga - log2 + 7) + A^o(a|Xi - X 2 |), (106) 

from which we see that f3\ = /3 2 and the eigenvalue condition is just 

- ci - v = -(log a - log2 + 7) + XoHXa - X 2 |), (107) 



20 



so that 

9 

a _ 2 e Cl+ly -T+ K "o(o|Xi-X2|) £ _ _ e -7r/2g+ci-7+-ft'o(a|Xi-X2|) (108) 

This expression for fc agrees with that derived for a pair of spirals in [12] by other methods. 

We observe that as the spiral separation |Xj — X 2 | — > oo, a approaches that corre- 
sponding to a single spiral of unitary winding number ()65p . 

Thus far we have matched at leading order and determined the eigenvalue a and 
corresponding asymptotic wavenumber k. To determine a law of motion of the individual 
spirals we need to match at first order. 



4.3.4 Outer limit of the first-order inner 

We sum the g-expansion of the outer limit of the first-order inner solution in the same 
way that we did for the outer limit of the leading-order inner solution. We first write 
equation (l8"oT) - ([8"7j) in terms of the outer variable to give 

-e^ ■ V/ = 6 2 V 2 /i - e 2 fi\V X o\ 2 ~ 2e 2 f V X o ■ V*i + fi~ 3fSfi, 
-/ie/ 2 ^ • V Xo = e 2 V • (/ 2 V X i) + e 2 V ■ (2/o/lVxo) + 2g/ /i - 4g/ 3 /i. 
We now expand in powers of e as 

Xi ~ +Xii(9,WH , 



to give 



h ~ /io(g,At) + e/u(g,/i) 



/io = 0, /u = -V£oo • Vxio, 



Motivated by the transformation we applied to £oo we write 

Xio = = , 

g/io ? 

whence (11091) becomes 

-^•V/? = V%, (110) 

where we have written xoo i n terms of ho- Recalling that ho = e qne ^ H (R) , and denoting 
as before 
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the left-hand side of (111 OR gives 

-f W ■ [ ^ e, + H (R)e q *e R j 

* ( -Vu 1 ^- - VuH^R)) - /ie^cos0 ( qni V 2e ^S- + V ie H' (R) 



fie** ev sin (p yqn e . s 

^e qn ^qn l s : m(P(R iqni - l A l e- iqril {yu - %Vu) + R' 1 ^ 1 ' 1 B^iVu + iV n )) 
- i2e qn ^qn e cos(j)(R iqne - 1 A e e- iqne (V 2e + iV u ) - R'^' 1 B^^ - iV u )) 

+ [tei*** ^ 1 ~ e 1 ) R- iqn t- 1 qn i B i e i ' ini (V u + iV-a) 



\iqnte qni<i> 



(e^R iqne A e e- iqne {V 2e + iV u ) - e^R^ iqne B e e iqnc {V 2e - iV u )). 



R 

Therefore, writing 

h = -nqniAte-^tyu + iV u )G 1 {R) e^+W - M n t B t ^{V 2t - iV u )G 2 {R) e^"^, 
gives 

G[ (gn e + ifd ■ x 

Gl + li + R? = R ' 

G' 2 (qn e - i) 2 G 2 n-iqni-i 
G * + li + R? = R 

with the general solution 



G x = - + 7ii? 1 -^ + y 3 R~ 1+iqne , 

Aiqni 

o—iqn^+l 

G 2 = + l2 R 1+iqn ? + 74 i?- 1 -^, 

Atqne 



where 71, 72, 73 and 74 are arbitrary constants. By comparison with the inner solution we 
find that 73 = 74 = 0; 71 and 72, will be determined by matching to the inner limit of the 
outer solution. Thus (with a redefinition of 71 and 72) the outer limit of the first-order 
inner solution is given by 

~ nA e e~* qn t(V u -iV 2e ) Rlcinf+1 e[qne+i)4> _ ^B^{V u + iV 2 i) R _ iqnf+l e{qm _ i)(f> 

1 4 4 

+ r )iR l ~ iqni e^ qne+i ^ + 72-R 1+il?r ^ e^ - *)^. (HI) 
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4.3.5 First order matching: determination of the law of motion 

To determine the law of motion we need to match the two-term inner expansion to the one- 
term outer expansion (in the notation of Van Dyke [T3], we impose (2ti)(lto) = (lto)(2ti)). 
Since we have summed the q expansion of the inner, but have only the first term in the 
q expansion of the outer, we again need to compare both series in terms of the outer 
variable. 

The outer limit of the two-term expansion in the inner region is 



X ~ - log h + ^ ~ - log (h + e/ii) . 

q qh q 



Thus, at leading order in q, we can perform the matching using h as in $31 Comparing 
(11 111) with ( 1971) we see that we need 71 and 72 to be of order one, in which case (with the 
usual expansion 7$ ~ 7^ + qja + • • • for 7 and V), and recalling (|93l)-(l94l) 



hi ~ 710 : -Re ^+ 720 : Re 



and in order to match with ( 1971) we require this to be equal to X • VG(X^) so that 

jlA eo e m ^ 2 (V uo -iV 2eo ) , G x (X e )-iG Y (X £ ) 



7io 



4 2 



flBo e e- in ^ 2 (V 1£ o + iV 2m ) G x (X e )+zG Y (Xi 
720 = 1 7. 



Now, writing (II 111) in terms of the inner variable r and expanding in powers of q we find 
that, as r — > 00, 

Xio ~ (^ o cos0 + V 2m sm<f)) + HfvG(X £ ) ■ e . (112) 
4.4 Solvability Condition 

We recall that equation ( 1851) is a linear equation of the form 

L(q,i/j )[^i] = b, 

where 

L(g,Vo)foM = (1 - ig)VVi + ^i(l -l^ol 2 )- (US) 
& = -^-V4. (114) 

In the case of nonzero g the operator L is no longer self adjoint. However, the adjoint 
operator is given by 

L*{qM[v] = (l+tq)V 2 v + v(l - |^o| V^ofoM* + ^) = L(-q,^o)[v\. (115) 



23 



Choosing v to be a non-trivial solution to the homogeneous equation L(—q, ipo) [v] = and 
using the Fredholm Alternative (integrating by parts) we obtain the solvability condition 



where D is a ball of radius r. The non-trivial solution of the adjoint equation are direc- 
tional derivatives Vt/'o ■ d of ipo, with q replaced by —q, where d is any vector in IR 2 . Since 
we are only interested in the leading-order law of motion, we expand the first-order inner 
equation ( jUSJ) in powers of q, and consider the solvability condition on the first term 
in this expansion. The leading-order (in q) operator is self-adjoint, and the non-trivial 
solutions are derivatives of the leading-order (e = 0, q = 0) inner problem. At leading 
order in q the solvability condition is exactly that of §2.2.11 namely 

^ V d) + * = 0. (117) 

4.5 Law of motion 

Finally, to obtain the law of motion we substitute (11121) into (11171) to give 

dX.p. 2nf ^_]_ 



dT Pep, 

where G is given by (1981) . 



rVGiXA (118) 



4.5.1 Law of motion for two spirals 

As an example consider two spirals at positions (A 1; 0) and (X 2 ,0), with Xi < X 2 . In 
that case 

VG(X 1 ,0) = (-P 2 aK' (a(X 2 -X 1 )),0). 
Using ji = 1/| loge|, q = vr/(2| loge|), and pi = p 2 the law of motion reads 

^ = (0, mi aK' (a(X 3 - A x ))) . (119) 

The direction in which each of a pair of spirals moves depends only on its own degree and 
not on the degree of the other spiral. This is exactly the opposite to what happens in fluid 
dynamics vortices, which would move depending only on the vorticity of the surrounding 
vortices. 

If we now take the limit q — > we would expect the law of motion (11191) to tend to 
the law of motion (1251) for q = 0. Since we scaled our outer region with q by setting 
q = 7r/(2| logeQ, letting q -> is equivalent to letting X 2 - X x -> 0. With X 2 - X t < 1 
equation (JUS]) becomes 

'0,-— ^— . (120) 



dT V X 2 -Xx 

We see that this does not agree with ( 125]) ; in fact, the two velocities are orthogonal to 
each other. 
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In order to match the law for spirals at the canonical separation with that for q = 
we need to consider an intermediate region, in which the separation of spirals is much less 
than the outer core radius. We analyse such a configuration in Sj5j We will find that the 
law of motion for such separations includes components both along and perpendicular 
to the line of centres, and matches with both the laws (j!19p and f[2"5l) in the appropriate 
limits. 



5 Interaction of spirals in the near-field 

We now assume that the spirals are separated by distances smaller than the canonical 
separation. The first step is to determine how this affects the eigenvalue a. If we expand 
( 11081) as |Xi — X 2 | — > we find a ~ |Xi — X 2 |~ 1 ^ 2 from which it seems that a should be 
larger than order one in the present scaling. However, we need to remember that we are 
redefining the e which appears in the definition of a. Since a is inversely proportional to e, 
this redefinition dominates the square root growth above, so that a is in fact algebraically 
small in e and exponentially small in q. 



5.1 Outer Region 

Proceeding as in §4.11 and expanding in powers of e as 

/ ~ fo(q,fj) + e 2 /i(<?,^) + ■ ■ ■ 

x ~ xo(q, + • ■ • , 

we find 

fo = 1, /i = -^|Vxo| 2 , (121) 

VXot = V 2 Xo + g|Vxo| 2 . (122) 

As before we write fi = jxq and assume (which will be justified a posteriori) that jx is of 
order one. Expanding 

Xoo . . 
Xo ~ 1" Xoi + qX02 H , 

q 

we find 

= V 2 Xoo + |Vxoo| 2 , (123) 
fiXoor = V 2 Xoi + 2Vxoo • Vxoi, (124) 

From and ijH we expect 

Xoo = C X {T) (125) 
to be constant in space, while the topological singularities are contained in xoi-, given by 

Xoi = C 3 (T) + Lrf, + C 2j (T) \ogR 3 + ^^R]) , (126) 

where C 2 j and C3 are time- dependent functions that will be determined by matching 
with the inner solution. As before, we use the notation Rj and <pj to represent the polar 
coordinates centred at the j-th spiral. Note that, since we are now working with x rather 
than h, it is clear that the corresponding wavefunction ip is single valued. 
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5.2 Inner Region 

The expansion in the inner region proceeds exactly as in §4.21 

5.3 Asymptotic matching 

5.3.1 Outer limit of the leading-order inner solution 

We again use the technique of Sections [3] and H] which allows us to sum all the terms in 
q in the outer limit of the inner expansion. The outer limit of the leading-order inner is 
again given by 

fo ~ 1 - -^|Vxoo| 2 H , 

Xo ~ Xoo(<?) H , 

with £oo = n e <j) + (1/q) log H (R) where, 

H = A e {q)e- igne R igne + B £ (q)e iqn ' R~ iqn \ (127) 

where 

Ai(q) ~ -Aio + An + qA i2 H , 

Q 

B e (q) ~ -B eo + Bn + qB l2 H . 

Q 

As before, comparing with the inner expansion gives 

A eo - B eo =0, (128) 
^1? In^ * = ~ W *>Tk = l,...,N. (129) 

A(q + D t o 

The remaining equations determining An and Bg will be found from matching with the 
outer region. 

5.3.2 Outer limit of the first-order inner 

The outer limit of the first-order inner solution is identical to that of §4.3.41 namely 

fi ~ -eVxoo • Vxio H , 



Xi 



where 



with 



hie- qxo ° 
Xw = 

q 



1 4 4 
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5.4 Inner limit of the outer 

To compute the inner limit of the outer we rewrite solutions (1125R and (11261) in terms 
of the inner variable X = + ex or equivalently Ri = |X^ — X| = er, and expand in 
powers of e. This gives 

Xo ~ — + C , 3 + n £ + C , 2 ,log(er) + G(X,) + eVG(X,)-x + O(e 2 ), (131) 

q 

where 

N c , 

G(X) = Mi + C 2j log |X - X,| + ^|X - X/, (132) 

a* 

We can now match this inner limit of the outer with the outer limit of the inner, given 
by (]127p . Since we have the full logarithmic expansion in the inner region, but only some 
terms in the logarithmic expansion of the outer region, we must write both expansions 
in terms of the outer variable before comparing terms. 

Expanding (11271) in powers of q using the relation = gives 

\og(A l0 (e- iqntloge + e iqnelost )/q) 
Xoo h n £ (f) 

q 

Ane~ iqnt log e + Bne iq1lt log e e ~^mi log e _ giqng log e 

+ int — : — ■, : — ■, — logi? (133) 



Aft) [e~ iqnt log e + e iqni log e ) e -iqn t log e _|_ e iqn e log e 

5.4.1 Leading order matching (lti)(lto) = (lto)(lti). 

We can now match the leading order terms by comparing (11331) with the leading order 
terms in e in the expression for the outer (11311) . We see that the leading term in ( 11331) 
is 1/glogl/g, which is not present in f 1 1 3 1 [) . but which is a constant (and constants are 
unimportant in x)- Thus we see that the outer expansion should really proceed aa3 

1 c 

Xo ~ -— + — + C 3 + n e <p + C 2i log(er) + G(X t ) + eVG{X e ) • x + • • • . (134) 
q q 

Then, matching H134[) and (11331) gives 

A eo {e- iqnelog€ + e iqnelog€ ) = e Cl , (135) 

A £ie -iqn e log e + Q^iqm log e 



= C 3 + G(X,), (136) 

e -iqn e loge _ e iqn e \oge 
ins : : : : = Co£- (137) 



A£ Q (e~ iqneloge + e iqneloge 

e -iqn e loge _ e iqn e \o; 
Q—iqrif log e _|_ ^iqrn log e 



1 this expansion can also be motivated by the single spiral solution {38] 
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5.4.2 First order matching (2ti)(lto) = (lto)(2ti). 

We now match the 0(e) terms in (I13ip with the 0(e) terms in the outer limit of the 
leading-order inner, which are 



Xio 



H A e e- iqnelose {V u - iV 2e )R 1+iqne e^ 
Aq ( y A£R iqne e~ iqnelogt + B(R~ iqni e iqniXog,L ) 

fi B e e iqnelost (V u + iV 2 e)R 1 ~ iqnt e- i * 
Aq ( y A£R i i ne e~ iqnelogt + BiR~ iqnt e iqntlo & t ) 
-f 1 R 1 - iqn e i<t> + 7 2j R 1+ign e-^ 



+ 



q ( y AeR iqril e~ iqniXo & t + BiR~ iqne e iqnelog( 



Expanding as q —> 0, 



,'A ll0 -flA m e- iqn ^(V li0 -iV 2£0 )\ t 

Yin ~ I -flC 

4 (A eo e- iqn * lo § e + B m e iqn i lo s e ) J 

±l2v-^B m e iqn ^(y UQ + iV 2 , 



4 (A m e- iqn ^ ogt + B eo e iqnelog( 
This must match with X ■ VG(X^), giving 



2J01 ) Re - i4> 



7io = ^ 1 2 [Ame + i%e 6 j , 

720 = ^ 1 2 + ^ J • 

Finally we are ready to write Xw back in terms of the inner variable r and expand in q 
to give 

Xio ~ — j- (^o cos + V 2 fl) sin 0) + — (y ia cos(</> - 2gn^ log e) + V 2eo sin(0 - 2gr^ log e)) 



-r 



cos(qn e loge){Gx0^e) cos(0 - gn^loge) + Gy(X £ ) sin(0 - gn £ loge)). (138) 



5.5 Law of motion 

The solvability condition on the first-order inner proceeds exactly as in $31 To obtain the 
law of motion we substitute (11381) into (j!17p to find 

dX oi 2 cos(gmloge) ± 

di /i sm(gri£ log ej 

Differentiating (I132f) gives 

AT / s 

v«(x) = E(,^ + C ^). 

Using (11371) this becomes 

sm(qrij\ log e|) e, 



dX, 


2 cos(gri£ 


log 


^1) 




rs_/ 

fx sm(qni\ 


log 


■A) 



V n ^ ^ (140) 



|Xj — X^| cos(gr?,j| loge|) |X 3 - — X^| 
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Since rij = ±1 for all j, we can simplyfy (11401) to 



dT 



11 iv^/ cos(gllogel) e r \ 

2n/g loge > n,,- ^ — 2_JZ 3 - - + - - , (141) 

1 ^ V J sin(g|loge|)|X,-X,| \Xj-Xt\J' 1 ! 



where we have taken /i = 1/| loge| as usual. We see that the law of motion in the near- field 
region interpolates between motion along the line of centres and motion perpendicular to 
the line of centres as q\ loge| varies from to 7r/2. 

5.5.1 Matching with the other limits 

If we take the limit as q\ loge| — > in (11411) we find 

, v N 

in agreement with (1251) . 

If we take the limit as q \ loge| — > 7r/2 in (11411) we find 

This should match with the limit v — > — oo of (11181) . As z/ — > — oo we see from (11081) that 
a — > 0, so that 

VG(X) ~ - X] 



|X-X 

Now, letting v — > — oo, a — > in ( 11 04ft gives, at leading order, 

TV 

i=i 

with solution a = e v / N , (3j = (3 = constant, for all j, in agreement with (jl08p . Thus in 
the limit as v — > — oo (11181) gives 



dX-z 



N 

dT -™< T, fx,-x,r 

j=i ,jjte 1 3 1 

in agreement with (11421) . 

6 Conclusions 

The main contribution in this paper is the description of complicated patterns of the 
complex Ginzburg-Landau equation with many moving spirals in terms of simple sets of 
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ordinary differential equations that provide a law of motion for the centres of the spirals. 
We have focussed on the case of spirals with winding number ±1 (one-armed spirals). 
Our results rely on the structural stability of such spirals; this is discussed in [7] where 
it is conjectured that 1-armed spirals are the most stable. 

Our analysis is based on the limit q — > 0, where q is a measure of the imaginary 
component of the coefficients in the equation. We find that the cases q = and q > 
behave very differently. For q > each spiral has an inner core (in which the magnitude 
of the wavefunction varies from zero to one), and an outer core (in which the level lines 
of the phase of the wavefunction vary from radial to azimuthal). The radius of the outer 
core can be related to the asymptotic wavenumber (the wavenumber at infinity) for a 
single spiral, and is exponentially large in q, tending to infinity as q — * 0. We find that 
the law of motion of spirals depends on the relative sizes of the separation and the outer 
core radius. Furthermore, the whole pattern oscillates at a frequency which varies slowly 
as the spirals move. Our main results can be summarised as follows. 



6.1 Laws of motion 

For < q <C 1, given a set of ±l-armed spirals with mutual separation of order 1/e, the 
spirals evolve on a time-scale given by T = e 2 t/\ loge| and satisfy the following laws of 
motion: 

(i) In the so-called canonical separation of the spirals, which corresponds to 
gj log e| = 7r/2 + qv + 0(q 2 ), with v = 0(1), the spirals' centres satisfy the law 
of motion 

«j ~ _^£, Wa | X< - Xj |)e,,, (143) 
where c\ = 0.098..., along with the set of linear equations 

N 

= ft( Cl + v - log a + log2 - 7 ) +'£l3 j Ko(a\X t - X,]), (144) 

for the parameters (3j, whose solvability condition determines a. This solvability 
condition shows that the corresponding asymptotic wavenumber, k = ea/q, is ex- 
ponentially small in q and evolves in time along with the positions of the spirals. 

(ii) In the so-called near-field separation of the spirals, which corresponds to inter- 
spirals separations such that < q\ loge| < n/2, the spirals' centres satisfy the law 
of motion 

N 



dT 



dX e cos(gloge) e r \ 

^ \ J sm (q loge X^ - X* X^-X^y 



(iii) In the limit as q — > the law of motion ( 11451) agrees with the one obtained by Neu 
in [9], 
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We can use the laws of motion in (i) and (ii) to form a composite expansion valid through- 
out the range < glog |e| < tx/2. We saw in §5.5. II that the law of motion in the overlap 
region (i.e. the limit of ( 11431) as v — > — oo, which is the limit of (1145ft as (/log |e| — > tx/2) 
is ( 11421) . Adding (1143j) and (I145p and subtracting (11421) gives the (additative) composite 
law of motion as 

dX £ Txn e a ^ Y 



o i i cos (g log e 

2n<glog e > w 3 - , 

sm(g loge 



sin(g| logeQ |X 3 - - X*| \Xj - X t 



N 

Alternatively, we can form the (multiplicative) composite expansion 

N 



EM X f -X 3 |)(,,|«e, +% ) (148) 



dX £ 2niq\og \e\a 

rsj 

dT dm 

which has the correct asymptotic limit in each region. These formulae are illustrated in 
Figure [U The additative composite is accurate in the region it was constructed to hold 
(g log | e| < vr/2), but the multiplicative composite remains a good approximation even 
for q log |e| > tx/2. 

We see that for small separations the interaction is along the line of centres, with like 
spirals repelling, opposites attracting. However, as the separation increases the direction 
of the interaction gradually changes, until it is perpendicular to the line of centres at 
large distances. 

While the attraction/repulsion of spirals depends on the winding numbers of both, the 
direction of rotation of one spiral about another depends only on its own winding number: 
the centre of positive spirals rotates around any other wpiral in an anti-clockwise direc- 
tion, while the centre of a negative spiral rotates in a clockwise direction about another 
spiral. Thus like positive spirals rotate in an anti-clockwise direction while separating, 
like negative spirals rotate in an anti-clockwise direction while separating, and unlike 
spirals translate while approaching. 

Since the motion is perpendicular to the line of centres in the canonical scaling, the 
question arises as to the existence of bound states. To answer this question the first- 
order correction to the radial velocity is needed. Our calculations indicate that the radial 
velocity remains of one sign, so that bound states are not possible for small q, in agreement 
with [2], who found that bound states are only possible for q > 0.845. 

6.2 Comparison with direct numerical simulations 

To compare with numerical simulations we write (11481) in terms of the original variables 
x and t to give 

g^q 2ri£qa r w/-i ( cos(g|loge| 



dt Bin 
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where a = ea satisfies 



= pA Cl - — - \oga + log2 - 7 +J2/3 j K (a\x i - Xj |). (150) 

v Zq J m 

In Figures [2] and [3] we compare the predictions of the canonical, near-field and com- 
posite laws of motion with a direct simulation of using second-order accurate finite 
differences in a square domain of side 800. The figures show the separation and rotation 
velocity of a pair of ng = 1 spirals when they are separated by a distance 60. 

We see that the multiplicative composite expansion captures the qualitative behaviour 
very well, and provides a reasonable quantitative prediction. 
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A Extension to non-zero b 

For clarity throughout we have only considered the case 6 = 0. However, the extension 
to general values of b does not introduce any further conceptual difficulty, it merely 
complicates some of the algebra. We indicate here briefly the necessary steps in the 
calculation. 

The outer equation for nonzero b reads 

e 2 (l - ib)^ T = e 2 VV + (1 + iq)(l - |^| 2 )^ - (151) 

q 

which in terms of modulus / and phase x becomes 

+ bfxr) = 6 2 V 2 / - e 2 /|V X | 2 + (1 - f 2 )f, (152) 
^\Pxt - bffr) = e 2 V ■ (fV X ) + qf(l ~ f) ~ e 2 -/ 2 . (153) 

q 

Upon expanding \ and / in powers of e we find that the leading-order terms satisfy 

/o = l, 

UXor = V 2 Xo + <?|Vxo| 2 - a 2 /q- 

where the parameter b does not appear; thus the rest of the outer calculation is identical 
to that of SU 

For b 7^ 0, the inner equations close to the £-th vortex are given by 

e^efr " ■ V/ + ebf XT - bf^ • V X ), 

= V 2 /-/|Vx| 2 + (l-/ 2 )/ (154) 

e ^ b ^F ' /v/ " ebffT + efxT ~ fd ~W ' Vx) 

e 2 o 2 

= V-(/ 2 V x ) + g(i-/ 2 )/ • (155) 

q 

Expanding in powers of e we find that at leading order 

= V 2 ^o + (l + ^o(l-|^o| 2 ), (156) 
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while at first order 

- Ml - ib)^ • Wo = V 2 ^ + (1 + zg)(^(l - |^o| 2 ) - MWi + WO), (157) 
or equivalently, 

-/i^-(V/o + fe/oVxo) = V 2 A-/ 1 |Vxo| 2 -2/ Vxo-Vxi + /i- 3/ 2 / 1 , (158) 

-/i^ • (/ 2 V X o - 6/oV/o) = V-(/ 2 Vxi) + V- (2/oAV X o) + 2qf f 1 -^f 1 . (159) 

Writing these in terms of the outer variable to find the outer limit of the first order inner, 

and ex panding in e as Xi ~ Xio(o)/* + Xn(<?) H and /i ~ f 10 (q) + e/n(g) H as in 

§4.3.41 we find the equation for xw is now given by 

- /x^ • Vxoo = V 2 £io - 2qf n = V 2 Xio + 2gVxoo • Vfto, (160) 

which again does not involve b. Thus the matching between inner and outer solutions is 
the same as in $H 

Finally we must derive the law of motion when b is not zero. To do so we start by 
defining the same inner product as before and consider the first order linear operator 
given in (j!57p . which is of the form 

L{q, ipo) [ipi] = w(ip , q, fi, b, dX e / dT) 
where L does not depend on b but where the non-homogeneous term is now 

Since the linear operator L is the same that that in §H the only place where b has an effect 
is in the solvability condition itself, which is obtained through the Fredholm Alternative. 
The new solvability condition is 

- /„« {"f 1 - ft ># ■ d ° = I s {f 1 - *> - 1&) } *■ < i6i > 

Since the left-hand side is O(fi), a nonzero value of b does not affect the velocity law at 
leading order, but it will alter the 0(/S) correction terms. 
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Figure 1: The velocity of rotation of a pair of +1 spirals as a function of q for e = 1/20 
and |Xi — X 2 | = 2 using the formulae for the canonical scaling (solid), the near-field 
scaling (dashed), the additative composite (dotted) and the multiplicative composite 
(dot-dashed). We see that the multiplicative composite is likely to be a much better 
approximation, especially for q > 7r/2 log |e|. 
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Figure 2: Rotation velocity of a pair of spirals 60 units apart as a function of q, using 
(11431) [solid], (11451) [dashed] and the composite (11481) [dotted]. The points correspond to 
a numerical simulation of ([I]). 
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Figure 3: Separation velocity of a pair of spirals 60 units apart using (11451) [dashed] and 
the composite (I148j) [dotted]. 
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